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We present the results of study of random sequential adsorption of linear segments (needles) on 
sites of a square lattice. We show that the percolation threshold is a nonmonotonic function of the 
length of the adsorbed needle, showing a minimum for a certain length of the needles, while the 
jamming threshold decreases to a constant with a power law. The ratio of the two thresholds is 
also nonmonotonic and it remains constant only in a restricted range of the needles length. We 
determine the values of the correlation length exponent for percolation, jamming and their ratio. 



PACS numbers: 05.40. +1, 64.60 Ak 
I. INTRODUCTION 

The problem of percolation is an old one but still 
new results appear and some unsolved questions re- 
main In general site percolation is defined on a d- 
dimensional lattice where each site can be either occupied 
with the probability c or empty with the probability 1 — c. 
Neighboring occupied sites form a cluster. If it is so large 
that it reaches the two opposite edges of the lattice, e.g. 
top and bottom, the cluster is said to be percolating. The 
lowest concentration of occupied sites for which there is 
a percolating (or spanning) cluster for an infinite lattice 
is called the percolation threshold c p . 

Another realization of the percolation problem is 
random sequential adsorption (RSA), in which objects 
(point particles, segments, rectangles etc.) are put on 
randomly chosen sites and the objects do not move 
It is also possible to consider RSA in a continuum j|. 

Jamming is a problem related to RSA percolation || . 
Again objects are placed randomly on the lattice sites 
until a concentration Cj is reached, where there is no 
room on the lattice for the next object. For point like 
particles Cj = 1, but for spatially extended entities c 3 < 1. 
Continuum models of jamming also exist 

The RSA models irreversible dissociation || and bind- 
ing of large ligands to polymer chains || . Another area 
of applicability is the deposition of large molecules on 
solid surfaces, like proteins JtJ or macromolecules on bi- 
ological membranes Q . The isotropic-nematic transition 
in the hard rods like polymers has been studied first by 
Flory |^] and later e.g. in jl0|]. Spatial organization of 
needles into a well organized nematic phase is however 
a different problem, not considered here. General forms 
of percolation models have a wide range of applications - 
from chemisorption, spatially disordered systems, porous 
materials, car parking and ecology ||, to separating the 
good and bad people at the entrance to Hades pT| . For 
overview of percolation, jamming and related problems 
see@. 

In a recently published paper |1J] Vanderwalle et al. 
studied the relation between the two transitions - per- 



colation and jamming. They used two kinds of objects 
- linear segments of length 2 to 10 and square blocks. 
They have found that the ratio of the two threshold con- 
centrations c p and Cj is constant c p /cj w 0.62, regardless 
of the length of the needle. 

In the present paper we extend the study of Vander- 
walle et al. to larger lattices and longer objects (we con- 
sider only linear segments). In particular we shall check 
the claim that the c p /cj ratio does not depend on the 
length of the needles. 



II. THE MODEL 

We consider a square lattice of size Lx L. On the sites 
of the lattice we put randomly linear segments (needles) 
of a given length a, with the constraint that the nee- 
dles cannot cross each other, although they may touch 
themselves. We used hard boundary conditions, i.e., the 
needles may touch the edge of the lattice but they can- 
not stick out of it - each needle must lay totally inside 
the lattice. Adopting open boundary conditions does not 
affect the results. 

To achieve simulation efficiency, our algorithm of de- 
position needles consists of two parts designed for two 
different regimes. Firstly when the current concentra- 
tion of the needles is small, we chose randomly, from a 
uniform distribution, the orientation (vertical or horizon- 
tal) and position of the upper left end of the needle to 
be inserted. If there is enough space on the lattice, the 
needle is deposited, if not we pass to the next try. After 
a certain number of adsorption trials we switch to the 
other regime where the dense routine is applied. A list 
of all empty sites and orientations still available is made. 
From that list a site is randomly chosen. We determine 
the direction of the needle and check whether the needle 
can be put there. In any case the site is removed from 
the list. The process is continued till the last item on 
the list. Such organization saves time, since we avoid to 
insert needles into densely packed regions. 

A cluster is defined as a group of sites linked by the 
needles. If there is an uninterrupted path between the 
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top and the bottom of the lattice, the cluster is said to 
be percolating or spanning, and the concentration of oc- 
cupied sites defines the percolation threshold c p . The 
concentration at which no more needles could be put on 
the lattice without violating the constraint determines 
the jamming threshold Cj. 

We have considered lattices of sizes L = 30, 100, 300, 
1000, 2500 and needles of length a = 1 .. 2000. On the 
smallest lattices only smaller needles were located. Av- 
eraging was done over 100 independent runs. We have 
checked that averaging over 1000 runs did not reduce the 
error (mean standard deviation - a) in a marked way. 



III. RESULTS 

Our main results of the simulations are presented in 
Figure 1, where the percolation and jamming thresholds, 
as well as their ratio, are plotted against the length of the 
needles (a = 1 .. 45). These data are obtained for lattice 
size L = 2500. As convergence and error analysis shows 
(see below) wc can safely accept them as the asymptotic 
(L — > oo ) values. 




FIG. 1. Thresholds for percolation c p , jamming Cj and 
their ratio c p /cj versus needles' length, a. Lattice size L = 
2500. Averaged over 100 samples. 

The percolation threshold for a = 1 .. 13 diminishes, 
then it begins to grow linearly with the slope 0.00071. 
The minimum value c pm i n = 0.463 is reached for a = 13. 
As seen in Figure 2, the a increases with the size of the 
needles starting from 0.001 (a = 1) up to 0.008 (a = 45). 
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FIG. 2. Percolation threshold c p versus needlel a \ength a. 
L = 2500, 100 runs. Short needles a = 1 ..45. 
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FIG. 3. Percolation threshold c p versus needles length a. 
L = 2500, 100 runs. Long needles a = 1 .. 2000. 

The increase of the percolation threshold for longer a 
is however quite clear. The appearance of this novel and 
unexpected feature is connected with the condition that 
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the needles may touch themselves but they cannot cross. 
In the simulations where the restriction has been lifted 
we observed no minimum but a monotonic decrease. In 
the model considered here the needles have the tendency 
to align in parallel not only with respect to the edges of 
the lattice but also to themselves (sec Figures 4 and 5), 
hence the needles form compact clusters. 




(a) 



FIG. 4. Snapshot of a spatial distribution of needles at the 
percolation threshold for L = 100, a = 5. 




(b) 



FIG. 5. Snapshot of a spatial distribution of needles at the 
percolation threshold for L = 100, a — 20. 

In case of horizontally oriented needles in order to move 
e.g. two steps down, two needles of length a are needed. 
The longer are the needles the higher is the percentage 
of occupied sites necessary for passing these two steps. 
The increase of c p (a) is to a certain degree compensated 
by vertically oriented needles, which however also form 
clusters, thus offering many equivalent ways for percola- 
tion. Further simulations for much longer needles indi- 
cate continuous increase in c p , although at a slower rate 
- see Figure 3. 

The jamming thresholds obtained from the simulations 
have much smaller error than that for percolation and 
even for a = 45 it is below 0.002. Values of cj, as a 
function of a, decrease according to a power law (very 
good fit for all a > 5) approaching the asymptotic value 



c* = 0.66 ± 0.01 (see Figure 6): 



Cj = c* + 0.44 • a-"- 77 . 
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FIG. 6. Jamming threshold Cj versus needles' length a on 
a log-log plot. L = 2500, 100 runs, a = 5 .. 45. 

The uncertainty of the exponent derived from the 
graph analysis equals 0.02. Clearly this behaviour dif- 
fers essentially from bare power law postulated in jjj : 

Cj ~a-°- 2 , (2) 

for the continuum (off - lattice) case of RSA of randomly 
oriented and highly anisotropic (length to wide) rectan- 
gles. Their a coincides with our length of needles a. In 
the discrete case we did not observe the maximum of Cj 
at a = 2 reported in M. The reason is that on the lat- 
tice the number of possible orientations of the needles is 
restricted to z/2 (where z is the coordination number of 
the lattice) in contrast to the continuum case. It is in- 
teresting that the asymptotic concentration for jamming 
(for a — > oo) is off lattice and it remains finite in the 
discrete case. 

Other interesting quantity in our model is the ratio 
Cp/cj as a function of a (see Figure 7). 



FIG. 7. Percolation to jamming thresholds ratio c p /c 



sus needles' length a. L 
a = 15. .45. 



2500, 100 runs. Logarithmic fit for 



It grows for a — 1 .. 3, then it stabilizes till a = 7 
and then it grows again. The plateau value of c p /cj w 
0.62 , the constant found in ]l2[ ). The growth for longer 
needles (15 < a < 45) could be fitted by a logarithmic 
dependence 



Cp/cj -0.50 + 0.13 -log a. 



(3) 
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FIG. 8. Percolation to jamming thresholds ratio c p /cj ver- 
sus needles' length a. L = 2500, 100 runs, a = 1 .. 2000. 
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Further simulations for longer needles (see Figure 8) 
support our claim of monotonic increase in c p /cj over 
wide range of a (even up to a — 2000) . We may conclude 
therefore that the universality claimed in [Q holds only 
in a rather restricted range of a € [3, 7]. As a matter of 
fact, the value of c p /cj for a > 7 shown in Table 1 of fl2| | 
is greater than those for a < 7 but the authors attribute 
it to the finite size effects. This is however most probably 
just the beginning of the growth of c p /cj. 

We analyzed the dependence of the obtained thresholds 
on the lattice size L and needles' length a focusing on 
convergence. It appeared that for the ratio a/L < 1/3 
the values of c p and Cj do not vary much with increasing 
L (keeping a constant) - see Figures 9 and 10. 
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FIG. 9. Convergence analysis of percolation, q^' jamming 
thresholds, Cj, and their ratio c p /cj, versus lattice size L. 100 
runs, a = 5. 



FIG. 10. 
thresholds, 
runs, a = 45. 



Convergence analysis of percolation, 

} j amming 

Cj, and their ratio Cp/cj, versus lattice size L. 100 



The error bars (here a) however decrease rapidly with 
L, whilst the difference of the thresholds for different 
lattice sizes is much smaller than the appropriate error. 
Thus it is safe to take the values of the thresholds from 
the simulations with L = 2500 as the asymptotic (exact) 
ones. 

The finite size effects can clearly be seen in Figure 3, 
where c p is drawn against a — 1 .. 2000 for L = 2500. At 
a = L/2 we can notice sharp change in the slope of the 
function c p (a). 

Consider now the dependence of a of c p , Cj , c p /cj on the 
lattice size, a is analogous to the quantity A defining in 
]l2[ the sharpness of the transition (non-percolating to 
percolating or non-jammed to jammed). Here however 
the power law approach to the asymptotic value p(oo) — 
p(L) - L~ x l v (ci. formula (3) in p| ) does not hold. 
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FIG. 11. Deviation analysis, a versus lattice(^ize L for 
several values of the needles' length. Percolation. 
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FIG. 13. Deviation analysis, a versus lattice^ize L for 
several values of the needles' length. Percolation to jamming 
ratio. 
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FIG. 12. Deviation analysis, a versus latticeC^ize L for 
several values of the needles' length. Jamming. 



We have found (see Figures 11 - 13) that the a for 
percolation (A p ), jamming (A,,-) and the c p /cj ratio (A r ) 
decrease with the lattice size according to the power laws 



lL -l/u r 



l/v p = 0.75 ±0.05, 
l/vj = 1.00 ± 0.05, (4) 
A r -L" i/l/ % l/v T = 0.77 ±0.05. 

Here v corresponds to the correlation length exponent j^j 

eHc-cpi-". (5) 

These values are, within the error bars, the same for all 
a = 1 .. 45 and agree with those found by Vanderwalle et 
al (if . Also Nakamura @ found Uj = 1.0 ±0.1 for RSA 



of square blocks. It seems therefore that the exponents 
v are good candidates for universal quantities. 
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(a) 



FIG. 14. Snapshot of a spatial distribution of needles at 
the jamming threshold for L — 100, a — 5. 



and 15 (jamming). Analysis based on examination of dif- 
ferent runs shows some regularity in the needles distribu- 
tion - we have found that the needles near the edges have 
the tendency to stick along the borders. Longer needles, 
for obvious reasons, form clusters of parallel alignment, 
as was already observed in fll2|. 



IV. CONCLUSIONS 

We have performed extensive simulations of RSA us- 
ing linear segments of size a = 1 .. 45 on square lattice 
sites. We have found that the percolation threshold is 
a nonmonotonic function of a, having a minimum due 
to parallel orientation of the needles, at a = 13, while 
the jamming threshold decreases to a non-zero constant 
with a as a power law. The ratio of the two thresholds 
is nonmonotonic too - after initial growth it stabilizes 
for some values of a, and then it grows logarithmically. 
Whether the asymptotic value is equal to one or below 
it is an interesting question. To answer it unequivocally 
is unfortunately beyond our computing power. The val- 
ues of the correlation length exponent v, for percolation, 
jamming thresholds and the ratio of the two, do not de- 
pend on the length of the needles and they are, within 
the error bars, equal to those found elsewhere [H plj 
for deposition of needles, rectangles or squares. 




(b) 



FIG. 15. Snapshot of a spatial distribution of needles at 
the jamming threshold for L — 100, a = 20. 

Examples of spatial arrangements of shorter (a = 5) 
and longer (a — 20) needles on a lattice 100 x 100 are 
shown in Figures 4 and 5 (percolation) and in Figures 14 
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